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In this letter we report our results in investigating edge 
effects of open antiferromagnetic Heisenberg spin chains with 
spin magnitudes S = 1/2,1,3/2,2 using the density-matrix 
renormalization group (DMRG) method initiated by White. 
For integer spin chains, we find that edge states with spin 
magnitude Sedge = S/2 exist, in agreement with Valence- 
Bond-Solid model picture. For half-integer spin chains, we 
find that no edge states exist for S — 1/2 spin chain, but 
edge state exists in S = 3/2 spin chain with Sedge = 1/2, in 
agreement with previous conjecture by Ng. Strong finite size 
effects associated with spin dimmerization in half-integer spin 
chains will also be discussed. 

PACS Numbers: 75.10.-b, 75.10.Jm, 75.40.Mg 



The antiferromagnetic Heisenberg spin chains has been 
a subject of immense interest in the last decade since 
Haldane conjectured that the low energy physics of in- 
teger and half integer spin chains are fundamentally dif- 
ferent [0. It is now generally believed that half-integer 
Heisenberg spin chains have gapless excitation spectrum, 
whereas gaps exist in integer spin chains (Haldane gap). 
More recently, there has been increasing interests in stud- 
ies of spin chains with defects |^-||. In particular, the 
properties of broken 5 = 1 quantum spin chains have 
received much attention because of the experimental ob- 
servation of 5 = 1/2 excitations localized at the ends of 
broken 5 — 1 spin chains More generally, one may 
address the question of whether edge states are genuine 
properties of finite quantum spin chains as in Fractional 
Quantum Hall Effect [^[(7). Recently, a theory of edge 
states based on the Non-linear-sigma model (NLcrM) 
plus topological 0-term has been developed by one of 
us js] where it was conjectured that edge states are gen- 
uine properties of antiferromagnetic quantum spin chains 
with spin value 5 > 1/2. In this letter, we shall ad- 
dress this question for both integer and half integer spin 
chains numerically using the recently developed density 
matrix renormalization group (DMRG) method IJ. We 
shall present results for open spin chains with spin values 
5=1/2, 1, 3/2, 2 up to chain length of 100 sites. 

The DMRG method has proved to be tremendously 
successful in studying 5 = 1 and 5 = 1/2 antiferro- 



magnetic Heisenberg spin chains H[[l| . The method was 
found to be particularly suitable for studying spin chains 
with open boundary condition and is thus well suited 
for our purpose of studying edge states. We use the in- 
finite chain algorithm jjj] in our study. Two new sites 
are added to the spin chain in each optimal step of the 
calculation from length L = 4 to L = 100. That is, 
open spin chains with even number of sites are studied. 
The number of the kept optimized states in our calcu- 
lation is m = 120. The largest truncation errors for 
S = 1/2,1,3/2, and 2 chains are found to be smaller 



than 10" 



io- 



3 x 10 and 10 respectively when 



the chain length reaches L = 100 in the final step. Prop- 
erties of the ground state and a few lowest excited states 
are obtained by looking at the lowest energy state with 
fixed total z-component of spin angular momentum 5' ot . 
In particular, the ground state corresponds to the lowest 
energy state in the sector 5' ot = 0. We shall look at the 
excitation energies of various states with different 5* ot 
and the corresponding average z-component of angular 
momentum on each site i =< S\ >. The dimmerization 
parameter q(i) =< 5j.5i_|_i — Si-i-Si > for the ground 
state will also be examined. 

We start with the excitation energies for integer spin 
chains. Fig.l shows the excitation energies E n — Eq for 
n = 1 to 3 as a function of chain length L for the 5 = 1 
spin chain. E n is the energy of the lowest energy state 
in the sector 5*°* = n. Energy is measured in units of 
Heisenberg coupling J. According to the valence bond 
picture Jl0| , ^2[ , for 5=1 spin chain, two 5=1/2 spins 
are left at two ends of the spin chain, and are coupled 
with effective coupling J e ff ~ Je^ L ^, where £ is the 
correlation length. For even spin chains, the coupling is 
antiferromagnetic and the resulting ground state is a spin 
singlet(5^ ot = 0). The lowest energy state with 5^ ot = 1 
can be constructed by exciting the singlet formed by the 
two edge spins into a triplet, and the excitation energy is 
of order J e f j , which goes to zero exponentially as length 
of spin chain increases. Excited states with larger 5'°' 
cannot be constructed by exciting only the edge spins any 
more and bulk excitations must be involved in construct- 
ing states with 5' ot > 1, implying that the excitation en- 
ergy will be of order E n — Eq ~ (n — 1) x Eh as L — > oo, 
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where Eh is the Haldane gap. This predicted behav- 
ior is clearly confirmed by our numerical result presented 
in Fig.l. The Haldane gap Eh and correlation length 
£ are estimated to be Ejj = 0.41J and £ = 6.0 from 
E 2 -E ^ 0.41J and E x - E ~ 0.7e~ L ^ ° for large L 
respectively, in agreement with result obtained by White 
Fig. 2 shows similar results for 5 = 2 spin chain. No- 
tice that in this case, the 5* ot = 2 state also has excita- 
tion energy going to zero as L — > oo, whereas higher 5* ot 
states have finite energy gap. This behavior is consistent 
with valence bond picture which predicts that for 5 = 2 
spin chain, the edge state spin magnitude is 5/2 = 1, and 
the two edge spins are coupled again by J e ff ~ Je~ L ^ y 
where £ is now the correlation length for 5 = 2 spin chain. 
In this case both 5* ot = 1 and 5' ot = 2 states can be con- 
structed by exciting edge spins only, and have excitation 
energy — > as L —> oo, which is exactly the behavior 
obtained in our numerical calculation. The Haldane gap 
and the correlation length are estimated to be Eh ~ 0.02 
and £ ~ 33 from E 3 -E ^ 0.02 J, E x -E ~ 0.1e- L / 33 , 
and E2 — Eo ~ 0.4e~ L / 33 for large L , respectively from 
our numerical results. However, the accuracy of these 
estimated numbers are much worse than the 5 = 1 case 
because of the much larger truncation error found in our 
calculation for the 5 = 2 spin chain. 

Next we turn to the excitation energies for half-integer 
spin chains. Fig. 3 shows the excitation energies E n — Eq 
for n = 1 to 3 as a function of inverse chain length 
L^ 1 for 5=1/2 spin chain. According to (abelian) 
bosonization theory |13| the excitation energy is given by 
E n —E = (5* ot ) 2 x (n/L). The bosonization prediction is 
also drawn on Fig. 3 (dash line). Notice that bosonization 
theory is correct only in the asymptotic limit L — > 00 and 
terms of order (LZn(L)) -1 are neglected It is thus 
not surprising that our numerical results and bosoniza- 
tion theory predictions show only qualitative agreement 
with each other. Notice also that bosonization theory for 
finite 5 = 1/2 spin chain predicts no edge states in this 
case. Fig. 4 shows similar results for 5 = 3/2 spin chain 
where the excitation energies are shown up to state with 
gtot _ 4 is apparent that E n — E ~ \/L in all four 
cases. Notice however, that edge state with spin magni- 
tude 5 = 1/2 is predicted to exist in this case, according 
to the conjecture by Ng Q. However it does not seem to 
show up clearly as for integer spin chains in the energy 
curve. 

To understand this behavior we examine the edge state 
in 5 = 3/2 spin chain in more detail. According to Ng 
H , the low energy physics of open 5 = 3/2 spin chain can 
be described by an effective 5=1/2 spin chain coupled 
antifcrromagnetically to two impurity spins with magni- 
tude Si m = 1 at two ends of the spin chain. The impu- 
rity spin will be partially screened by a Kondo type ef- 
fect at low energy ||, leaving "free" spins of magnitude 
Sedge = 1/2 at the ends, coupling ferromagnetically 
to the bulk 5=1/2 spin chain. Again, the coupling 
of the "impurity" spin to the bulk 5=1/2 spin chain 
can be analyzed using renormalization group technique. 



For ferromagnetic coupling, the corresponding operator 
is marginally irrelevant, implying that a free Sedge = 1/2 
is left at each end of the spin chain as L — > oo. For a 
finite spin chain, an RKKY type coupling between the 
two edge spins will be found, and the resulting ground 
state is a spin-singlet for even chains, as is in the case 
of integer spin chains. The RKKY coupling Jr between 
the two edge spins comes from exchange of spinwaves, 
and has a length dependence Jr ~ g/L, where g is some 
effective coupling constant. Thus we expect that the ex- 
citation energy for the edge spins E e d is proportional to 
g/L. Notice that at large distance L, the coupling con- 
stant g will be renormalized with g — > g / (1 + gln(L)) 
(ferromagnetic Kondo effect) and corresponding correc- 
tion to energy E e d will be found. In particular, at large 
L, E e d ~ 1 / (Lln(L)) . Assuming that the low energy bulk 
excitations are described by bosonization theory as an ef- 
fective 5=1/2 spin chain, we conclude that the lowest 
(5' ot = 1) excited state for 5 = 3/2 open spin chain is an 
edge excitation as in integer spin chains, since the bulk 
excitation energies scale as 1/L, and will be always higher 
in energy then E e d as L — + oo. More detailed analysis 
of the energy spectrum supports our edge state energy 
analysis which we shall discuss in a later paper. In the 
following, we shall present a more direct numerical evi- 
dence supporting our edge states picture for Heisenberg 
spin chains. 

Fig. 5 shows the expectation value < S l z > for the 
5*°* = 1 state for all four spin chains 5 = 1/2, 1, 3/2, 2 
with 100 sites, i.e. i = 1, 100. The four cases are ar- 
ranged from top to bottom in increasing order of spin 
magnitude, i.e. the top one is for 5 = 1/2, and bot- 
tom one for 5 = 2, etc. For spin singlet ground states, 
< SI > will be zero for all sites i. Thus < S\ > for the 
gtot _ 2 state gives information about the wavefunction 
of the first excited state in the spin system. Notice that 
for all four spin values, Fig. 5 shows that the excited states 
all carry staggered magnetization because of the under- 
lying antiferromagnetic interaction. For 5 = 1/2 spin 
chain, Fig. 5 suggests clearly that the first excited state 
can be thought of as a standing spin wave. The same 
qualitative result is also obtained from bosonization the- 
ory of finite 5=1/2 spin chain. In fact the feature can 
be understood in the much simpler 5=1/2 XY-modcl 
with J z = 0. The 5=1/2 XY- model can be mapped 
onto a non-interacting spinless fermion model. It is easy 
to show that similar feature exists in this case. The in- 
troduction of nonzero J z term just modified the feature 
quantitatively. It is clear from Fig. 5 that the wavefunc- 
tion of the first excited states in the 5 > 1/2 spin chains 
are qualitatively different from that of the 5=1/2 spin 
chain. In fact, in all three cases we considered, < 5* > 
suggests clearly that the first excited states are all con- 
sist of spin excitations which are localized around edges 
of spin chains, i.e. they are edge excitations. It is also 
clear from the figure that the correlation length of the 
edge state in 5 = 2 spin chain is much larger than the 
corresponding correlation length of the 5=1 spin chain, 
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in agreement with general expectation and with results 
obtained from excitation energy analysis. One can also 
study < S\ > for the 5*°* — 2 states and looks at the 
difference in < S\ > between the 5f* = 2 and S\ ot = 1 
states. For the S = 3/2 spin chain, the difference is 
found to have a spinwave like feature similar to < S\ > 
for S = 1/2 spin chain with 5' ot = 1, suggesting clearly 
that the spin magnitude of the edge state is Sedge = 1/2, 
and excitations with 5* ot > 1 involves bulk excitations, 
in agreement with previous conjecture by Ng g. 

We have also examined the dimmerization parameter 
q(i) =< Si.Si+i — Si.Si-i > for the ground states of the 
four spin chains we considered. For integer spin chains, 
q(i)'s are found to be nonzero around ends of spin chains 
and decay exponentially as one move into the interior, in 
agreement with prediction by Ng For half-integer 

spin chains, <?(i)'s are found to be large at the ends and 
decay with power law into the interior of spin chain. Sim- 
ilar results were also obtained by White || for 5=1 and 
5=1/2 chains. For 5=1/2 chain, the dimmerization 
parameter q(i) can again be examined using bosonization 
techniques where it can be shown that q(L/2) ~ (1/L)", 
where /3 = 1/2 in bosonization theory. Similar behavior 
is observed in our numerical results for both 5=1/2 
and 5 = 3/2 chains, with ~ 1/2 for S = 1/2 chain and 
(3 ~ 0.6 for 5 = 3/2 chain. The huge difference in be- 
haviour of q(i) between integer and half-integer open spin 
chains reflects the fact that half-integer spin chains are 
more susceptible to spin-Peierls instability and furnish 
another interesting parameter distinguishing between in- 
teger and half- integer spin chains. 

Lastly, we want to list our calculated ground state en- 
ergies (per site) eo for the four spin chains. We obtain 
e = -0.443147 for 5=1/2 chain and e = -1.401484 
for 5=1 chain, in good agreement with exact Bethe 
ansatz solution and the result obtained by S.White et.al 
|TTf respectively. We have also computed the ground 
state energies for 5 = 3/2 and 5 = 2 chains with good 
precision, with eo = —2.8283 for 5 = 3/2 chain and 
e = -4.7608 for 5 = 2 chain. 

Summarizing, in this paper we have performed nu- 
merical studies on low energy properties of finite quan- 
tum spin chains with spin values 5 =1/2, 1,3/2, 2 using 
the Density-Matrix Renormalization Group techniques 
which has been applied successfully to study 5 = 1/2 
and 5 = 1 spin chains. We find that all the results we 
obtained are consistent with the edge states picture con- 
jectured by Ng m. In particular, we present for the first 
time numerical evidence for existence of edge states in 
5 = 3/2 and 5 = 2 finite quantum spin chains. The 
dimmerization parameter q(i)'s for the ground states are 
also studied where the huge difference in behavior be- 
tween integer and half- integer spin chains is pointed out. 
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FIG. 1. The excitation energy E n — Eo for even S = 1 
chain as function of chain length L. The black circles, black 
squares, and black triangles are for E\ — Eo, E2 — Eo , and 
E3 — Eo respectively. 

FIG. 2. The excitation energy E„ — Eo for even S = 2 
chain as function of chain length L. The open circles, black 
circles, black squares, and black triangles are for E\ — Eo, 
E2 — Eo, E- A — Eo, and E4 — Eo respectively. 

FIG. 3. The excitation energy E n — Eo for even 3 = 1/2 
chain as function of inverse chain length 1/L. The black 
circles, black squares, and black triangles are for E\ — Eo, 
E2 — Eo, and E3 — Eo respectively. The dashed lines are 
En - Eo = n 2 ty/L. 

FIG. 4. The excitation energy E n — Eo for even S = 3/2 
chain as function of inverse chain length 1/L. The open cir- 
cles, black circles, black squares, and black triangles are for 
Ex — Eo, E2 — -Eo, E3 — Eo, and E4 — Eo respectively. 
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FIG. 5. The configuration < S\ > for the lowest states of 
S = 1/2, 1,3/2,2 chains with S* ot = 1 and length L = 100. 
The zero point for S = 1/2,1,3/2,2 chains are shifted to 
4, 3, 2, 1 respectively. The circles, black circles, black squares, 
and black triangles are for the < SI > of chains with 
5 = 1/2, 1, 3/2, 2 respectively. 
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